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The heaviest members of the SM - the gauge bosons W, Z and the top quarks and antiquarks - 
may form collective bag-like excitations of the Higgs vacuum provided their number is large enough, 
both at zero and finite temperatures. Since Higgs vacuum expectation value (VEV) is significantly 
modified inside them, they are called "bags". In this work we argue that creation of such objects 
can explain certain numerical studies of cosmological baryogenesis. Using as an example a hybrid 
model, combining inflationary preheating with cold electroweak transition, we identify "spots of 
unbroken phase" found in numerical studies of this scenario with such W — Z bags. We argue 
that the baryon number violation should happen predominantly inside these objects, and show 
that the rates calculated in numerical simulations can be analytically explained using finite-size 
pure gauge sphaleron solutions, developed previously in the QCD context by Carter, Ostrovsky 
and Shuryak (COS). Furthermore, we point out significant presence of the top quarks/antiquarks 
in these bag (which were not included in those numerical studies). Although the basic sphaleron 
exponent remains unchanged by the top's presence, we find that tops help to stabilize them for a 
longer time. Another enhancement of the transition rate comes from the "recycling" of the tops 
in the topological transition. Inclusion of the fermions (tops) enhances the sphaleron rate by up 
to 2 orders of magnitude. We finally discuss the magnitude of the CP violation needed to explain 
the observed baryonic asymmetry of the Universe, and give arguments that the difference in the 
top-antitop population in the bag of the right magnitude can arise both from CP asymmetries in 
the top decays and in top propagation into the bags, due to Farrar-Shaposhnikov effect. 



I. INTRODUCTION 

A. The Baryogenesis 

The question how the observed baryonic asymmetry 
of the Universe was produced is among the most difficult 
open questions of physics and cosmology. Sakharov [T] 
had formulated three famous necessary conditions: the 
baryon number and the CP violation, with obligatory de- 
viations from thermal equilibrium. Although all of them 
are present in the Standard Model (SM) and standard 
Big Bang cosmology, the known mechanisms creating the 
baryon asymmetry can only produced effects many orders 
smaller than the observed amount, usually expressed as 
the ratio of the baryon density to that of the photons 
ns/n 7 ~ 10~ 10 . To find a scenario in which this puzzle 
can at least in principle be solved is one of the motivation 
of this paper. 

It follows from general anomaly relation and, in partic- 
ularly, from known electroweak instanton and sphaleron 
solutions, that baryon number can be violated by certain 
processes providing a change in topology of the gauge 
field. In the broken phase of the SM we live in today these 
processes are basically frozen, by a very high barrier as- 
sociated with Higgs VEV. On the other hand, in the sym- 
metric phase - the electroweak plasma - the sphaleron 
rates are suppressed only by the powers of the coupling: 
those rates are high enough to wipe out any primordial 



asymmetry, if it has been generated cosmologically prior 
to the electroweak transition era. These consideration 
force us to search for resolution of the baryon asymmetry 
puzzle at narrow temperature interval at or right below 
the electroweak phase transition. For review of the field 
see e.g. [2H3- 

The bubbles associated with the expected 1-st order 
electroweak phase transition had attracted significant at- 
tention in 1980's-1990's. Moving walls of such bubbles 
may supply the necessary local deviations from thermal 
equilibrium. The rate of thermal sphaleron transitions, 
both in equilibrium and near the bubble walls, has been 
studied extensively, see [SHE] and many more. However, 
as the experimental limits on the Higgs mass have over 
time evolved upward, to Mh > 116 GeV today, it became 
clear that the first order transition is actually impossible 
in the SM. (Transition is a crossover for Mh > 80GeU, 
see again [5] for details and the original references). Af- 
ter this has fact has been acknowledged, people either 
looked at various phenomena beyond the SM such as its 
supersymmetric extensions, which may still allow some 
window for the 1st order transitions. 

An alternative is the modification of the standard cos- 
mology. Instead of the usual Big Bang scenario, with its 
adiabatically slow crossing of the phase transition and 
tiny deviations from equilibrium, the so called "hybrid" 
scenario has been proposed [22l [23] , see also [2H [25] and 
many subsequent works. It combines the end of the 
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inflation era with the establishment of the electroweak 
broken phase. In it the nonperturbative production of 
long-wavelength gauge and Higgs fields happens far from 
equilibrium, returning to equilibrium later, at a temper- 
ature well below the critical temperatures T c of the elec- 
troweak transition; therefore it is sometimes called the 
"cold" electroweak scenario. While based on some fine 
tuning of the unknown physics of the inflation, it avoids 
many pitfalls of the standard cosmology, such as "era- 
sure" of asymmetries generated before the electroweak 
scale. 

In such a scenario there are coherent oscillations of the 
gauge/scalar fields [26l|27]: those have been studied in 
detail in real-time lattice simulations [2"9l [50] . Some of 
these authors focuses on the problem of generation of pri- 
mordial magnetic fields, see e.g. [28] in the same scenario, 
which we will not discuss. The simulated models include 
two scalars - the inflaton and the Higgs boson - and the 
electroweak gauge fields of the SM, in the approximation 
that the Weinberg angle is zero (Z is degenerate with 
W) . All fermions of the SM are ignored: the effect of the 
top quark in particular is the subject of the present pa- 
per. After inflation ends, all bosonic fields are engaged 
in damped oscillations for relatively short time, at the 
end of which the Higgs VEV and gauge fields stabilize 
to their equilibrium values, with the bulk temperature 
Tbuik ~ 50 Gey, well below the critical (crossover) tem- 
perature. 

In this work we will use the results of these simula- 
tions as a quantitative example, although some our re- 
sults should also be applicable for standard cosmology 
as well. We start with explaining some of the results of 
these simulations (the sphaleron size and rate) in simple 
analytical models, and then proceed to discuss what are 
changes induced by the top quarks. 

Another ingredient of the baryogenesis problem, the 
CP violation, had also evolved during the last decades. 
The nature of observed CP-odd effects, first discovered 
in kaon decays and lately in decays of the B mesons, 
is rather firmly established and its origins traced to the 
complex phase of the CKM matrix. It has been argued 
that it should enter as the so called Jarlskog determinant 
, with all mass differences explicitly present polynomially, 
leading to effect ~ 10~ 19 , way too small for baryogenesis. 
The magnitude of the CP violation in the SM is however 
still hotly disputed: the effect must be proportional to 
the Jarlskog determinant only if all fermionic masses ap- 
pear explicitly in the numerator. This would be the case 
if the relevant scale of all loops is much larger than all 
the masses. This is not the case e.g. in the famous CP 
violation in the kaon decays. It has been argued recently 
[55] that next-to-leading order dim-6 operator appears 
which contains 1/toJ? in the coefficient: if so, numerical 
study of the possible role of this operator [5H] found CP 
effect even 4 orders of magnitude larger than needed! We 
will return to critical discussion of this issue at the end 
of the paper. 

Vast literature exists on possible leptogenesis scenar- 



ios, generating the asymmetry in the leptonic sector, for 
overview see e.g. [5J. If this is the case however, leptonic 
asymmetry should still be converted into the baryonic 
asymmetry, again at or near the electroweak transition 
scale. We will not discuss this option in this paper, to- 
gether with many other ideas beyond the SM. 

B. The W - Z -top bags 

Generic argument for multi-particle states due to Higgs 
attraction is that however weak the forces can be, if it is 
an universal attraction, a sufficiently large number N > 
N c of particles it will become very strong. For example, 
an extremely weak gravity force binds the planets and 
stars. This generic argument however assumes that Higgs 
is light enough, so that sufficiently many particles can be 
collected inside the region of a size < 1/Mjj, so that the 
universal Higgs-induced forces are now yet cut off by its 
mass. In the context of Higgs attraction it was discussed 
in particularly for light quarks in [9], in which case the 
iV c is astronomically large and it worked only for very 
light Higgs bosons. 

More recent round of ideas has focused, not surpris- 
ingly, on the heaviest member of the Standard Model 
(SM), the top quarks. The first question was whether 
a single top quark would significantly modify the Higgs 
vacuum in their vicinity to make a "bag" [TUJ E] . This 
idea was soon refuted (see e.g. |12j ) by quantum one- 
loop corrections. Similarly, the idea does not work for 
"realistic" Higgs mass and binary it or tt systems. 

Frogatt, Nielsen and collaborators [T5J were the first 
to look at multi-top system. They provided simple 
hydrogen-atom-like estimates for a system of N = 12 top 
quarks (6 top and 6 anti-top for the lowest S±/2 orbital), 
which suggested that its binding energy can be large and 
nearly cancel the mass! Unfortunately, our more accurate 
mean field calculations [14 also refuted it. The binding 
is nonzero for massless Higgs, although it is much smaller 
than estimated in [13]. With the "realistic" Higgs mass 
we found that the 12 top-antitop system is unbound. Fur- 
ther discussion of the binding of the 12 top system can 
be found in |15j . who refined the conclusion in the form 
of the maximal Higgs mass at which binding occurs: 

Mh < M c (12) » 49 GeV (1) 

Systematic study of the top-bags in relativistic regime 
has been started in the paper [14]. One-loop quantum 
corrections to the bags have been investigated in a sep- 
arate paper |16j . following basically the line of thought 
of Farhi et al. We have found that at Yukawa coupling 
corresponding to the top quark the quantum (one-loop) 
effect are still under control, for any AT, while it is not 
so for hypothetical fermions few times heavier than tops. 
Significant progress has been made in the companion pa- 
per [17] . in which we have studied binding of various 
bags, filled with bosons (W,Z) and fermions (i,t). We 
have found that bosonic bags may exist even in vacuum, 
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for "realistic" Higgs masses ~ 100 GeV , while purely top 
bags do not. Those bosonic bags should however have 
very large number of quanta, of the order of thousands, 
to get bound. 

One reason to study the W - ^-top-bags is method- 
ical: they are only the third (and highly relativistic) 
family of manybody bound objects, besides atoms and 
nuclei. Unfortunately their experimental production at 
accelerators seems to be impossible, mostly because it is 
highly improbable to produce many tops in a small vicin- 
ity from each other. In this paper we suggest another 
motivation for their study, the cosmological one. One 
of the main messages of this paper is that metastable 
(mechanically stabilized) multi-quanta bags replace the 
"metastable bubbles" considered 20 years ago, being the 
arena in which the baryon number violating sphaleron 
transitions (and possibly also CP violation) takes place. 



C. Sphalerons with and without the Higgs 

Semiclassical description of the tunneling through the 
barrier, separating topologically distinct gauge fields, is 
given by the famous electroweak instanton solutions, 
which however leads to extremely low tunneling rate per 
unit time and volume, normalized by the only scale of 
momenta, the temperature T, to a dimcnsionless quan- 
tity 

TtunneUng/T 4 ~ exp(-4n/a w ) ~ 10-™ (2) 

At finite temperatures another option appears, the so 
called sphaleron transition, which is due to a thermal 
excitation onto the barrier. The original electroweak 
sphaleron solution found by Klinkhamer and Manton 
(KM) [5] has the energy Ekm (the height of the bar- 
rier) of about 14 TeV . So , at the temperatures at the 
electroweak transition scale T ~ T c ~ 0.1 TeV, the corre- 
sponding Boltzmann factor is prohibitive 

Tkm/T 4 ~ exp(-E KM /T) <~ lO" 60 (3) 

However, this estimate is too naive, as the parameters of 
the electroweak theory are strongly renormalized near T c . 
More accurate calculation of the equilibrium sphalerons 
rates at the electroweak cross over region give much larger 
rates. For recent update see e.g. [18] who estimated those 
in the range 

T/T 4 m 10~ 20 (4) 

including rather large preexponent calculated in Refs. [BJ- 
8 . Such rates are however still too small for the solution 
of the baryogenesis puzzle. 

The main reason of the increase of the rate, from the 
KM sphaleron to revised one is of course the reduction 
of the Higgs VEV from its vacuum value v = 246 GeV to 
near-zero. (Recall that there is no Higgs-related barrier 



in the symmetric phase, and thus no exponential sup- 
pression: but this is not a blessing since a transition rates 
faster than cosmological expansion leads to "wipe out".) 

In the hybrid scenario one asks: what is the barrier 
height and the sphaleron transition rates for a no-Higgs- 
VEV bag of a specific size pi Although with a completely 
different (QCD) motivation, the analytic answer to this 
question was already known, given by the so called COS 
sphaleron |19j . The main obstacle for finding it earlier 
was that, lacking a nonzero Higgs VEV, classical gauge 
theory is conformal and has no dimensional parameter. 
Thus the minimal sphaleron mass can only be defined if 
some additional condition is imposed, e.g. that the r.m.s. 
size of this object is fixed 

9 J B(ry <Pr V ' 

Adding another constraint, that the Chern-Simons num- 
ber of the configuration Ncs is also fixed to some value, 
one can minimize the energy over all possible magnetic 
field configurations and find a "sphaleron path" , the set 
of configurations which describe the barrier separating 
one topological valley from the next. 

COS result for the energy as a function of the Chern- 
Simons number can be written in a parametric form, with 
the parameter k = 

E C os= 3 -^(l-k 2 r (6) 
9 9 

Ncs =\sign{k){2 + \k\){l-\k\) 2 (7) 

The k = 0,Ncs = 1/2 configuration is the top of the 
barrier -the COS sphaleron. 

We will be arguing below that the topological objects 
found numerically in real-time are in fact close to COS 
sphalerons with the total energy 

Ecos ~ 2 TeV « E KM * lATeV (8) 

much smaller than in the broken phase. (Once again: 
it is only possible because they are located in the "no- 
Higgs spots" with depleted Higgs VEV.) Obviously the 
corresponding Boltzmann factor is dramatically reduced, 
falling into the range in which one can discuss baryogen- 
esis. 



D. The aims and the structure of this paper 

In a sentence, we will show that the presence of the 
W - Z-top bags should significantly enhance the rate of 
baryon number violating processes. 

A bit more detailed summary is that there are at least 
three important effects we will consider: 
(i) As we just discussed, the baryon number violation 
rates inside such bags are enhanced by many orders of 
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magnitude due to the absence of the Higgs VEV and 
related high barrier. In other words, instead of KM 
sphalerons one should use the COS ones, reducing the 
barrier height from 14 to only about 2 TeV. We will be 
able to analytically estimate the sphaleron rates, which 
compare well with numerical simulations. 

(ii) Accounting for top "bags" in this setting we find 
that tops help mechanically stabilize the bags. In a 
way, the metastable bags play a role similar to that of 
the bubbles (present if the electroweak transition be 
of the 1st order), namely to enhance deviations from 
equilibrium for rather long time. 

(iii) Last but not least, the "recycling" of top quarks, 
present in the bag, effectively lower the barrier further, by 
about 600 GeV. CP odd effects may lead to top-antitop 
population difference in the bags, which will result in 
asymmetric diffusion in the baryon number. 

Now about the structure of the paper. We start in sec- 
tion [Ii] with the discussion of the the results of numerical 
simulations |29| and then qualitatively explain them us- 
ing the COS sphaleron. After that, we turn to top and 
W bags in section |TTT| We discuss their production and 
lifetime. Finally, in section IV we will study how the 
sphaleron transition itself is changed, with the account 
for fermions. 



II. HYBRID INFLATION SCENARIO: BOSONIC 
SIMULATIONS AND THEIR DISCUSSION 
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(i) One important finding of the simulations is that 
the initial coherent oscillations of scalars soon give way 
to the usual broken phase. The most important feature 
is persistence of "no-Higgs spots" in which Higgs VEV is 
very far from the equilibrium value v and is instead close 
to zero. The gauge fields in them have however rather 
high magnitude. Figfl] (from [25]) show an example of 
a snapshot of the Higgs field modulus. Typically the 
volume fraction occupied by such "no-Higgs spots" is of 
the order of several percents and is decreasing with time. 

(ii) The second important findings is that of topologi- 
cally nontrivial fluctuations of the gauge fields. As shown 
in the lower Fig{T] those are well localized. According to 
the anomaly relation, changes in topology leads to vio- 
lation of B + L (baryon and lepton) numbers, via emis- 
sion/absorption of fermions, 12 in the SM. 

It was observed that topological fluctuations happen 
only inside the "no-Higgs spots" mentioned above. In- 
deed, this becomes apparent from the distribution of the 
topological charge shown in the lower part of Fig[T] for 
the same time configuration of the Higgs as shown in the 
upper part of the Fig[l] . Of course, it is only one snap- 
shot, but the authors found from the simulations that it 
is true for the whole sample. 

The fraction of no-Higgs spots which induced topolog- 
ical transitions is also in the range of a percent. More 



FIG. 1: (Color online). The contour plots from 29 of the 
modulus of the Higgs field |</>(a;)| 2 (upper plot) and the topo- 
logical charge density Q(x. t) at time mt = 19, for the model 
Al, 7V S = 48. Red (dark) areas in the upper plot correspond 
to small VEV, while yellow (light) bulk corresponds to the 
broken phase. On the lower plot lumps of the topological 
charge density appear as red regions (dark in black and white 
display). While most of the no-Higgs spots do not have the 
topological transitions, all transitions seem to be inside the 
spots. 



precise measure is the so called sphaleron rate is defined 
by the mean square deviation from zero of the Chern- 
Simons number 



r(t) 



i 



dAN* s 



W dt 



(9) 



(here and below all quantities are defined via one charac- 
teristic mass parameter m: for the simulations its value 
is about 264 GeV, close to v. For the definition see Ap- 
pendix A.) 

Since the process only exist for finite period of time, 
ref.[2H] reports its time integrated value defined by some 
time integral 
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during which it basically happen. (Too early there are 
no gauge fields and too late there are no no-Higgs spots.) 
By the end of the process (mt < 45) it is in the range 



7-10" 



(11) 



for more details about different parameters sets see Ta- 
ble II of [25]. This quantity, as well as of course snap- 
shots like those shown in Fig|l] directly give the spatial 
distance between the topological fluctuations fl^m = 
20..30. 




FIG. 2: (Color online). From [55]. Two snapshots of the 
topological charge, at times mt = 18 and 19. 

(iii) Space-time evolution of the topological charge Q 
is shown in Fig[2] as two snapshots, and in more detail 
in Fig(3]for the maximal B 2 (t),E 2 (t) and Q(t), for one 
particular fluctuation. One can see that the fluctuation 
at some time moment is very much concentrated in a 
small spherical cluster (the upper one in Figj2]) is followed 
by an expanding spherical shell (the lower one in Figj2]) 
which gets near-empty inside. 

The decomposition into electric and magnetic compo- 
nents of the field is shown in Fig[3] One can see that 
the fluctuations starts as nearly (90%) magnetic object 
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FIG. 3: (Color online). From [29]. Upper: The spatial profile 
of the charge density peak for various times. Lower: The value 
of the electric, magnetic and charge density at the center of 
the peak as a function of time. 



at mt ~ 17, with the electric field and thus the topological 
charge ~ EB peaking some time later. Then there ap- 
pears an expanding shell, followed by the magnetic field 
rebounding to its secondary maximum of smaller ampli- 
tude. We will return to discussion of all those features in 
the next chapter. 

Some interesting findings of those studies we will not 
discuss but yet would like to mention. In Ref.[31 it has 
been shown that "hot spots" have some topology in terms 
of the Higgs winding number, and that their favorite val- 
ues is about 1/2 (thus named half- knots). There is some 
correlation between this winding number and Ncs (and 
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thus the sphaleron formation) which deserves to be stud- 
ied further. 

Before modelling these numerical results, let us sum- 
marize what is it exactly we would like to achieve: 

(i) to model the structure of the "no-Higgs spots" by 
WZ-top bags 

(ii) to model the topological fluctuations in them by the 
COS sphalerons, explaining their structure, size and ul- 
timately the rate 

(iii) to discuss specifically the modifications coming from 
the top quarks, as those were not included in numerical 
simulations. 



A. "No-Higgs spots" as the WZ bags 

In Ref.|17j we have made some first steps toward un- 
derstanding of the states of the gauge bosons in cer- 
tain spherically symmetric Higgs backgrounds. Our main 
finding was that while outside of the bag the gauge 
quanta W, Z are massive and have 3 polarizations, in 
the bags the longitudinal component obtains additional 
strong repulsive potential. It also is mixed with electric 
polarization, which is indirectly affected by this potential. 
Only magnetic one remains unaffected by this repulsion: 
their properties were documented in j!7j . 

In principle, at finite T one needs a complete set of 
W, Z states for the partition function of the bag. Fur- 
thermore, in order to model the "spots" appearing near 
reheating, one should also do it out of equilibrium, which 
is even more complicated. However, since such states 
are not yet known in full, we will start here with two 
grossly simplified limits, assuming that those spots con- 
tain so large number of W that they can approximately 
be treated macroscopically. 

Thermodynamically stable bubbles of the symmetric 
phase can only coexist with the bulk in the broken phase 
only at the transition temperature T c , provided the phase 
transition is 1st order. As we already mentioned in the 
Introduction, such option has been excluded by the com- 
bination of experimental limits on Higgs mass and lattice 
studies of the transition. And we also consider a "cold 
model" in which the equilibrated bulk is at a tempera- 
ture well below the critical one, T < T c . Therefore, we 
discuss the objects which are not thermodynamically but 
only mechanically stable. 

What we consider here are a very simple variant of the 
bag model. Let us assume that rescattering are rapid 
enough to make kinetic equilibrium inside the bag. There 
is no chemical equilibrium, as the number of Ws pro- 
duced are determined by early time dynamics: so we 
need to introduce both the internal temperature Tj„ and 
the fugacity £w = exp(u w /T in . Rescatterings will force 
the internal momentum distribution to be thermal 



r d k 

N w = V spotgw J 



1 



(12) 



(2tt)3 ^exp(e(k)/T in ) - 1 
even if there is no chemical equilibrium, say the num- 



ber of gauge degrees of freedom can be different from 
the maximal number gw = 6 (e.g. only magnetic modes 
excited, without electric or longitudinal ones, as we dis- 
cussed in the companion paper) or £w different from 1. 
Condition of mechanical stability we will write as 



Pin = B + p bu 



Ik 



(13) 



with the bag constant B created by the Higgs potential, 
assuming zero VEV inside 



B = Xv 4 /A = mV/4 



(14) 



(loop and T-dependent corrections can be easily included 
but are ignored. All other degrees of freedom of the SM 
such as quarks and leptons are ignored here because they 
have not yet been produced.) Using further Boltzmann 
approximation (ignoring 1 in denominator of (12)) and 
also ignoring pbuik for the estimate, one gets the mechan- 
ical stability condition in a form 



B = 9w£,w~ 
or the internal temperature 

l/4/gWx-l/4 . 



n 2 T?„ 



45 



Mm « 174 GeV 



(15) 



(16) 



So, if it is an equilibrium W gas (all factors are 1) the 
internal temperature inside the spot is indeed well above 
the electroweak critical temperature T c ~ 100 GeV, mak- 
ing it a spot of symmetric phase. 

Since the W bag is only mechanically stabilized, it is 
relaxing by cooling and thus disappear relatively quickly. 
Indeed, this is what happens in simulations, giving the 
characteristic lifetime of the no-Higgs spots 



That spots ~ 20/m 



(17) 



It is however important for the sphaleron rate that - in 
the bag approximation we use - the inside temperature 
Ti n is not changing while this bag shrinkage takes place, 
because it is related with the Higgs bag constant. 



B. The sphalerons, their structure and size 

Trying to understand the results of these numerical 
simulations, one can ask in particular the following ques- 
tions: 

(i) How so high-magnitude gauge field can be produced? 
Why does it happen only inside the no-Higgs spots? 

(ii) What is the energy needed for topological fluctua- 
tion? What is the gauge field structure? Can one under- 
stand its subsequent evolution and decay products? 

(iii) Can one estimate the rate of their production? 
The fact that topological transitions happen only in- 
side the no-Higgs spots is rather simple to understand: 
as we already mentioned in the introduction in the bro- 
ken phase the height of the barrier (the mass of KM 
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sphaleron) is prohibitively large. In a spot the height 
is not zero, however, as it has finite size p. As we already 
mentioned, it is given by the mass of the COS sphaleron 
~ 1/p (see @). 

In a scale-invariant classical YM theory the size p is 
an arbitrary parameter. In a Big Bang however ( as 
well as in numerical models we discuss) this size is to 
be determined by the optimal scenario maximizing the 
transition rate. 

Although the bag production is a result of a compli- 
cated nonequilibrium process, it is reasonable to assume 
that large-size bags are cut off exponentially in bag vol- 
ume. Furthermore, we model it by a thermal fluctuation 
with the inside temperature. In a simple bag model the 
pressure is p = Pi n (T) - B while the energy density is 
e = €in{T) + B. Mechanical stability condition for the 
total pressure of the bag p = (more accurately, the out- 
side pressure in the r.h.s., which is very close to zero ), 
or Pin(T) = B. Furthermore, if the inside of the bag is 
filled by massless gas of any kind, e in (T) = 3pi n (T) = 3B 
and thus the energy of the bag is 4BV. Since our bags 
are not asymptotically large, perhaps some coefficient (to 
be later determined) can be used instead: so we finally 
write the bag Boltzman factor as exp(-j3BV) 

So, the optimal size is the extremum of the following 
expression 



Fsph(p) ~ exp 



Sir 2 



(18) 



Here the first term schematically represents the proba- 
bility to create large bags, while the second is the Boltz- 
mann factor for the COS sphaleron, with the mass from 
Q. The resulting optimal spot size is then 



U0T in g 2 B) 



1/4 



(19) 



independent on the degrees of freedom filling the bag. 
Comparing it with the "numerically observed" pm ~ 3.9, 
Figs . 2 and 3, at the moment of maximal B, one may 
extract the value of the parameter /? which turns out a 
bit smaller than 1/Ti„. We repeat that the model is crude 
and we only need it to discuss later the influence of the 
tops on the whole picture. 

It is useful to discuss at this point how semiclassical 
are those sphalerons. Note that the parametrically large 
quantity here is electroweak coupling 4ir/g 2 = l/a ew , 
which determines that the sphaleron consists of many 
gauge quanta and thus can be treated semiclassically. As 
the radius contains its power 1/4 we get only factor 2 
from it. The number of gauge bosons involved can be 
estimated from the action/ft which is about O(10). It is 
perhaps still large enough for the semiclassical analysis 
we use, but - as we will see later - not large enough to 
ignore the back reaction from the 12 fermions associated 
with the electroweak anomaly. 



C. The shape and field structure of the sphaleron 

Ideally the sphaleron transition may proceed with the 
total energy of the gauge quanta exactly equal to to the 
height of the barrier. In this case at the time the system 
reaches the maximum - the "sphaleron moment" - the 
kinetic (electric) energy is zero, after which the system 
may fall downhill into a topologically different configura- 
tion. 

An interesting option pointed out recently by Kuchicv 
[331 134] is that the optimum energy may be less than the 
height with tunneling through the barrier's top. That 
however cannot happen in classical simulations we dis- 
cuss. 

One may think there would be some extra energy 
needed to fall over the top in those simulations. And in- 
deed, as seen in Fig(3j this is the case. The time evolution 
of the magnetic and electric fields have their "sphaleron 
moment" , defined as the maximum of B 2 (t). One indeed 
finds that the kinetic-to-potential E 2 /B 2 energy ratio is 
indeed small at this moment, of the order of percents. 

Does the magnetic field fits well to COS sphaleron so- 
lution? The profile of the the magnetic field in COS 
configuration is given by the following simple expression 



g 2 (r 2 + p 2 y 



(20) 



which (unlike the KM one) is just spherically symmetric. 
This form does indeed fit well the observed shape of the 
B 2 at the sphaleron moment. The maximum of B 2 /g 2 
can be related with its radius, yielding mp « 3.9, which 
is very close to optimal size we got above. This value 
corresponds to the total energy of the COS sphaleron 



Etot = 3ir 2 /g 2 p « 2 TeV 



(21) 



As we already mentioned in the Introduction, it is 7 times 
less than the KM sphaleron mass, and for the tempera- 
tures we are dealing with T in = 200 - 100 GeV it makes a 
huge difference for the rate. 

The sphaleron Boltzmann factor can now be estimated, 
as we know both the total energy and the internal tem- 
perature T in from our bag model for the no-Higgs spot: 

exp(-E tot /T m ) « ea;p(-2000 GeV/174 GeV) « 10~ 5 

This is not yet the end of the estimate, since semiclassical 
sphaleron rate has also a significant preexponent. It has 
not been calculated for COS sphaleron yet, so we use the 
KM one 



Vm 4 2ttto 4 



N tr N rot ( 



which includes the unstable frequency oj_ w 2M\y , as well 
as the numbers due to translational modes N tr = 26 and 
rotational modes N rot » 5300. There are also factor ~ 
0(1) from the non-zero mode determinants. We used 
here a 3 = g 2 /4ir = ^2M w /g 2 w T. 



Combining all the factors we find that numerical value 
of preexponent and exponent nearly cancels out, leaving 
crudely 



r/Vm 4 ~ 10" 



(24) 



with accuracy say an order of magnitude or so. With 
that accuracy it agrees with the results of the simulations 
which also finds that the number of sphaleron transitions 
per spot is indeed about several percents. 

What should happen after the sphaleron moment? 
Sphaleron decay is classical rolling of the classical (high 
amplitude) gauge field downhill, from the (sphaleron) top 
into the next classical vacuum . This process was exten- 
sively studied numerically for KM sphaleron. Remark- 
ably, an analytic solution of the time-dependent explo- 
sion of COS sphaleron has also been found, see the COS 
paper. We will not describe the details of that here, just 
give the expression for the late time t » p profile of the 
energy density of the expanding shell. 



47re(r, t) 



8tt 



g 2 p 2 r 2 



(l + (r-t)2/p 2 ) 



(25) 



Comparing the explosion of COS sphaleron with nu- 
merical data one can see both the similarities and the 
differences between them. As seen from Figj3j there is 
an empty shell formation at some time. However the 
inside of the shell does not remain empty: in fact the 
topology and magnetic field have a secondary peak (of 
smaller magnitude). Qualitatively it is easy to see why 
it happens. The COS sphaleron is a solution exploding 
in zero Higgs background, with massless gauge fields at 
infinity. In the numerical simulations we discuss such ex- 
plosion happens inside the finite-size cavity. As the gauge 
bosons of the expanding shell hit the walls of the no-Higgs 
spot, they are massive outside. With some probability 
they get reflected by this wall back. This is the reason 
why a secondary splash (the second magnetic peak) is 
produced. It would be possible and interesting to study 
this problem separately, obtaining COS-like explosion in 
a spherically symmetric bag: we hope to do so elsewhere. 

Let us make a small theoretical digression here. For 
both for KM and COS sphaleron explosions there was 
well known controversy about relation between the 
Chern-Simons number and the baryon charge for non- 
static solutions. At time going to infinity the Chern- 
Simons number during the sphaleron explosion has been 
found to be different from naively expected ANcs = 1/2. 
This happens because even a weak field at late time still 
preserves some topology. Nevertheless, it was somewhat 
troublesome before the issue of the fermion number viola- 
tion has been resolved by Shuryak and Zahed in [32] , who 
found the analytic solution of the Dirac eqn in a time- 
dependent exploding COS background. They found that 
any fermions starting from the sphaleron zero mode are 
indeed excited into physical (positive energy) states at 
late time. The outgoing fermions have simple momen- 
tum spectrum 



which is close to thermal with T e f t = 2/p, and the average 
quark energy < E q >= 3/ p. Thus the amount of baryon 
number violation by sphaleron transition is in fact as 
expected from the anomaly. 



III. THE ROLE OF THE TOP QUARKS 

A. Top quark production and concentration in the 

bags 

Since the original numerical simulations have included 
the gauge fields but ignored fermions, we have to discuss 
first, at quite qualitative level, what their effects can be, 
relative to that of the gauge fields. Those tops would be 
added to the metastable bubbles of the symmetric phase, 
the no-Higgs spots, like the W discussed above. However, 
there are at least three important differences: 

(i) their fugacity is expected to be larger, due to more 
rapid production rate of t relative to W 

(ii) the degeneracy factor gj t = 12 is twice larger than for 
gauge bosons (6). 

(iii) their binding to Higgs bags is several times larger 
than for W, which makes it easier to satisfy the mechan- 
ical stabilization of the bag. (For details on that see [17] .) 

The first question is at what time and how the top 
quarks would be produced. Top quark has the largest 
coupling to Higgs: so it is produced first via HH -* it 
process. Let us crudely compare its rate relative to that 
of HH -* WW 



T(HH tt) 
Y(HH -> WW) 



m w 



20 



(27) 



where in the first estimate follows from the fact that the 
(vacuum) masses are proportional to the corresponding 
coupling constants, which appears in the power 4 in the 
rate. This simple estimate suggests that top production 
may be by about an order of magnitude more rapid than 
that of Ws, and thus they would in fact dominate in the 
formation of "spots" . 

Obviously, specific simulations are needed to test top 
production rates. We are not currently able to eval- 
uate top production, give out-of-equilibrium and from 
strongly fluctuating Higgs fields at the preheating stage. 
As a crude guess one may look at an equilibrium density 
in the bulk before the appearance of broken phase VEV 



c(tt) 



4tt 2 



T A « LIT 6 



(28) 



n L (k)=p(2kp) 



2 -2k p 



(26) 



The next step is concentration of top quarks into the 
bags. The simple reason this is happening is that they 
are massive in the broken phase and massless in the sym- 
metric one. As discussed in detail in [T7], their binding 
can be as large as 100 GeV or so. Less trivial reason is 
that a "Id kink" of the Higgs field (a 2-d surface where 
Higgs VEV crosses zero) possesses fermion zero modes. 
In practice it means that top quarks can glue themselves 
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to the = surfaces and be transported along them, 
eventually into the remaining island of the symmetric 
phase. The binding energy in such case is nearly all the 
top mass, 172 GeV, several time the bulk temperature. 
It indicates that the effect can be rather robust: but only 
dedicated simulations can decide how effective this mech- 
anism can actually be, and what fraction of the produced 
tops end up in the bags. 

In order to see possible degree of such concentration, 
note that after tops are collected into the bags, they re- 
main effectively massless their, so their density inside is 
given by the same expression, with high Tj„ ~ 200 GeV 
needed to balance the bag pressure. The density in the 
bulk is given by the integral with the top mass and low 
T bulk ~50GeV 



n bulk ~ 9tt 



JM 



dEE\J E 2 - M 2 
2^2 



exp(-E/T) (29) 



which is smaller than the density n max inside the spots 
by about factor 300! It is thus clear, that after the broken 
phase is established in the bulk, most (~ 99%) tops should 
either be annihilated or collected into the bags. 

Let us first estimate it from above, assuming that all 
heavy quarks can be collected into the remaining spots 
of the symmetric phase. If so, the number of tops per 
spot is just the ratio of their densities 



jv(tt) = !*&L ~ io 2 (-) 3 

'"spots *"* 



(30) 



where T is a characteristic temperature of the Higgs field 
at early time, T/m ~ 1. Thus we expect N t+ i of the 
order of a hundreds of tops+antitops be collected into a 
no-Higgs spot. 

Since we are not able to model the process quantita- 
tively, in a highly fluctuating Higgs background we will 
introduce a parameter, normalizing this process to that 
of the W. So, we will denote the number of tops normal- 
ized to Ws which made into the "spots" : 



jyspot 

NT 



(31) 



The population of the spots including W, tops and anti- 
tops would then be ( 1 + 2k) times larger than in bosonic 
simulations. 



B. The W and the top quark lifetime in the bag 

The lifetime of the W according to their weak decays 
is, in our units, 



m 

Twm ~ 50 

Yw 



(32) 



The lifetime for bound W's is presumably reduced, ac- 
cording to their total energy increasing by another factor 
2 or so. This timescale is several times larger than the 



lifetime of the "bags" in the numerical simulations, which 
is of the order of 10-20 in the same units. 

Thus the weak decays can be neglected (as they were 
in the numerical simulations) and the observed lifetime 
in the simulations should be attributed to evaporation of 
the W's into the bulk, by thermal excitation. Indeed, the 
"hot spots" of the reduced Higgs VEV were considered 
there as mere fluctuations of the fields. We however look 
at them as metastable objects, which can come into me- 
chanical equilibrium and decay more slowly due to heat 
transport as well as strong and weak interactions, ignored 
so far. 

As discussed in [T7] in detail, we found that tops tend 
to be on the surface of the bag, due to the influence of 
the "zero mode" phenomenon. For how long would i, t 
quarks remain in the bag, before their decay? Because 
of hierarchy of the couplings, one should first look at an- 
nihilation via strong interactions. One of such processes 
leads to production of other quarks via it — > qq, another 
is the annihilation into gluons it — ► gg, which give the 
lifetime tj 4 ~ 200/T in , see Appendix B. This is about an 
order of magnitude longer time than the lifetime of the 
bosonic "no-Higgs spots" in [29] , 

(Another process it — ► HH has a bit smaller coupling 
and less final states, especially in color, so it is clearly 
subleading.) 

After strong annihilation of pairs, the remaining num- 
ber of tops or antitops are due to random charge fluctu- 
ations during the formation stage 



Ntiv 



(33) 



Assuming that this is the case, one may ask how the 
tops would decay further. In vacuum the physical top 
quark is much heavier than gauge quanta, so they decay 
t -*■ b + W + in the first order in weak interactions. The 
corresponding width is well known 



m, 



16 



( * )(l + 2 — 



(34) 



' w 



However, for both t and W bound in the bag it is by 
no means certain whether nit > mw or not. As we had 
shown in [T7] , for large enough bags the lowest top levels 
are in fact lower than the lowest W ones. 

When the corresponding top energy levels are still 
above the W states, such weak reaction is possible. Crude 
estimate of the weak decay width can be obtained if in 
this formula we substitute masses by energies, writing 
(for E t > Ew) the width as 



Et 



io-%i -E w /E t y 



(35) 



In practice, this leads to a lifetime of the order of T wea k ~ 
m/T~ 10 3 . 

When the top quarks are heavier than Ws in the bag 
decays t -*■ Wb become impossible. The weak decays 
would the proceed into the 3 body, say t -> blvi like in 
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the usual beta decays. This produce an extra power of 
a w and much stronger dependence on the energy released 
because of 3 body phase space. This width, crudely esti- 
mated as 

^~M^-J (36) 

where by AE t we mean the change in the total energy 
when one top disappears, is way too small and exceed all 
the scales in the problem (except the time of cosmological 
expansion) . 

In summary, mr 4 ~ 200 is the lifetime due to to strong 
top annihilation, while their weak decays take no less 
than mt ~ 10 3 - 10 time, or even much more if the bags 
are large enough. This should be compared to Twm of 
about 100 for the W's, which would clearly be shorter 
and thus limit the lifetime of the bags. The lifetime rm 
of about 20 for the bosonic spots in simulations is due to 
thermal evaporation: presumably a deeper binding and 
decrease of internal temperature due to tops will improve 
it. We conclude with a crude estimate of the bag lifetime 
as rm ~ 100. 



the exponent small, within the level of uncertainties of 
the model used. 

Now we turn to the discussion of the preexponent. Ob- 
viously it depends nontrivially on the presence of the tops 
in the bag, as fermions are related to the probability of 
topological fluctuations by the anomaly: we will turn to 
this effect in the next section. If the T 4 in the rate comes 
from some thermal average, one may think that it will 
come multiplied by (at least) one power of DOF factor. 
If so, again the increase of the DOF in the bag due to 
tops would be partially canceled out by a reduced inter- 
nal temperature, since the product is the bag pressure B. 
The preexponent simply has to be recalculated, for COS 
sphalerons in the bag environment (in which many fields 
such as gluons are absent because of different cosmolog- 
ical scenario). 

We have argued above that tops have longer lifetime in 
the bag, which may increase the overall bag lifetime. At 
the other hand, bag stability relies on W, Z presence and 
thus it cannot exist longer than the W lifetime, r^m ~ 
200. Thus the potential enhancement of the hot spot 
lifetime cannot be larger than factor 10. 



IV. THE SPHALERON TRANSITIONS IN THE 
TOP-STABILIZED BAGS 

A. The exponential effects and the lifetime 

As argued before, tops would be an important in- 
gredient in the metastable bags with near-zero Higgs 
VEVs. Now we turn to the next question: how much 
the sphaleron rate can be affected by their presence, as 
compared to purely bosonic ones in the simulations 29J? 

Since in this case there are more particles in the bags, 
with the tops added to gauge bosons, one may naively 
think that the bag energy would grow. But using a simple 
bag model we have already pointed out that this would 
not be the case: the mechanical stability require total 
pressure to be zero, which related the internal pressure 
to the bag constant B. No matter which light degrees of 
freedom are inside, the bag energy would scale as BV. 
The Boltzmann factors for sphalerons however would still 
depend on the number of degrees of freedom Ndo f in the 
bag because the temperature depends on it. This again 
follows from the mechanical stability condition 



N, 



DOF 



T7L ~ B 



(37) 



Thus adding tops to the bag simply reduces the internal 
temperature. Returning to the exponential approxima- 
tion for the sphaleron rate (181, we now put into it the 



optimal radius ( 19 ) and also rescale the temperature as 
in (37). One finds that the exponent of the sphaleron 



rate depends on the number of degrees of freedom as 

~ (Ndof) ' 25 and changes by just few percents only. 
We thus conclude that the cancellations due to mechan- 



ical stability condition (37) makes the effect of tops on 



B. Recycling the top quarks 

The well known Adler-Bell-Jackiw anomaly require 
that a change in gauge field topology by AQ ± 1 must 
be accompanied by a corresponding change in baryon 
and lepton numbers, B and L. More specifically, such 
topologically nontrivial fluctuation can thus be viewed 
as a "t'Hooft operator" with 12 fcrmionic legs. Particu- 
lar fermions depend on orientation of the gauge fields in 
the electroweak SU(2): since we are interested in utiliza- 
tion of top quarks, we will assume it to be "up" . In such 
case the produced set contains t r tbtgC r CbC g u r UbUgT (i, e , 
where r, b, g are quark colors. We will refer below to this 
process as the -»■ 12 reaction. Of course, in matter 
with a nonzero fermion density many more reactions of 
the type n -*■ (12 - n) are allowed, with n (anti) fermions 
captured from the initial state. 

Although in this work we try to follow thermal lan- 
guage (local kinetic equilibrium) as much as possible, let 
still discuss what should be done in a purely dynami- 
cal out-of-equilibrium setting. In principle, one should 
proceed quantum mechanically, starting and ending with 
certain number of ingoing and outgoing quark and gauge 
quanta in certain in and out states (e.g. fixed momenta) 
projecting those into the fermionic and gauge field con- 
figurations of the semiclassical theory of the sphaleron 
process. The sphaleron itself is just one - saddle - point 
at t = on the path, which start from convergent waves 
and end with divergent ones. 

The whole classical solution describing the expansion 
stage at t > has been worked out for COS sphaleron ex- 
plosion, and for the "compression stage" at t < one can 
use the same solution with a time reversed. At very early 
time or very late times t -*■ ±oo the classical field become 
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FIG. 4: Schematic picture of fermionic level motion as a func- 
tion of time. Open and closed circles indicate unoccupied and 
occupied states. 



weak and describe convergent/divergent spherical waves, 
which are nothing but certain number of colliding gauge 
bosons. Fermions of the theory should also be treated 
accordingly. Large semiclassical parameter - sphaleron 
energy over temperature - parametrically leads to the as- 
sumption that total bosonic energy is much larger than 
that of the fermions, so one usually ignores backreac- 
tion and consider Dirac eqn for fermions in a given gauge 
background. For KM sphaleron and effective T we dis- 
cuss, this parameter would be ~ 70, which is indeed large 
compared to 12 fermions. However in the case of COS 
sphaleron we are going to use the number is about ~ 10, 
comparable to the number of fermions produced. It im- 
plies that backreaction from fermions to bosons is very 
important. 

To our knowledge, the only (analytic) solution to Dirac 
eqn of the "expansion stage" was obtained in [32] , it de- 
scribes motion from the sphaleron zero mode at t = and 
all the way to large t -*■ +oo, the outgoing physical state 
of fermions, with momentum distribution (26). A new 



element we are adding now is that its time-reflection can 
also describe the compression stage, from free fermions 
captured by a convergent spherical wave of gauge field at 
t -*■ —oo and ending at the sphaleron zero mode at t = 0. 

How the corresponding fermionic levels move is 
schematically shown in Fig[4] In the Dirac sea notations, 
the level moves from some negative initial energy at early 
time to positive one at late time. For COS sphalerons 
the mean energy is E q = -3/ p - to E q = +3/ p By symme- 
try, the level is crossing zero at t = 0: the corresponding 
fermionic zero energy state is known as the sphaleron zero 
mode. If the level was occupied (as most of levels below 
the surface of the Dirac sea are) it remains occupied and 
one gets new fermions produced. If there was antiquark 
(an unoccupied state or a hole, indicated by the open cir- 
cle) it will remain unoccupied, thus having no meaning 
after the level moves to positive energy, where most levels 
are unoccupied. So, antiquarks - holes at negative energy 
sea - are actually decelerated by the radial electric field 
to zero energy modes, their energy being transmitted to 
the gauge field. 

This conclusion is very important: it implies that en- 



ergy of the initial fermions can be incorporated and used 
in the sphaleron transition. If initially fermions were in 
certain state (e.g. certain level of the Higgs bag \in> ), 
the fermion part of the sphaleron amplitude would in- 
clude the amplitude of their capture into a "sphaleron 
doorway" state \door way > by a convergent gauge field 
wave. As we will discuss in detail below, we will find that 
the optimum way to generate sphaleron transition is to 
use 3 initial top quarks, considering the 3^9 fermion 
process instead of the original -* 12 one. In order to 
satisfy Pauli principle and fit into the same sphaleron 
zero mode, colors of the 3 quarks of each flavor should 
all be different. 

Let us show by simple estimates that under conditions 
we are discussing the -*■ 12 fermion production pro- 
cess is actually impossible. The final kinetic energy for 
fermion produced by COS sphaleron is E q = 3/p which is 
« 200 GeV for the mp = 3.9 example displayed as typi- 
cal in Figj3] Multiplied by 12 fermionic species produced 
it would require 2ATeV of energy which apparently ex- 
ceeds the total available energy of the gauge fields in the 
topological fluctuations observed numerically. 

(Note finally, that inclusion of the CP-violating pro- 
cesses - such as induced by CKM matrix - may lead to 
a small difference between recycling of the top and the 
antitop quark amplitudes. This asymmetry would then 
be a mechanism for the baryogenesis.) 

The 3^-9 fermion process saves a lot of energy, as in it 
the initial top quark energy can be completely transferred 
from the "sphaleron doorway state" to the gauge field 
during the compression stage. In the example we discuss, 
with mp = 3.9, this mean energy of 3 top quarks is 



AE = 3* 3/p* 600 GeV 



(38) 



An estimate with exponential accuracy we find its en- 
hancement by the factor of about 



* recycling 



exp(A£/T) ~ 20 



(39) 



using the hot-spot temperature Tj n . 

Note that the same amount of energy should also not 
be produced in the final state: thus in terms of energy 
balance this amount is doubled, saving ~ 1200 GeV and 
making the transition possible. Yet this gain does not 
increase the rate any further, as the system falls downhill 
classically the probability of everything to happen being 
just 1, due to unitarity. 

Can one continue to use the same mechanism for more 
initial-state fermions, for example considering the cap- 
ture of 3 tops plus n other quarks, gaining another factor 
~ e nAE / T ? In particular, should one focus on 6 -*■ 6 tran- 
sition, which is "energy neutral"? The answer to this 
question is negative: for quarks other than t there is no 
reason to get trapped inside the no-Higgs spots, as Higgs 
effect on them (the mass) is too small compared to even 
low bulk Tout and they occupy the bulk of the volume. 
Roughly their density is thus smaller than that of top 
quarks by the factor (T out /Ti n ) 3 , which is few percents. 
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The cube of this factor obviously cannot compensate the 
gain from the exponent. 

The prccxponent for the 3 9 fermion process should 
include projection of the initial states of fermions (say 
bound states of ttt) in a bag onto the appropriate 
"fcrmionic doorway state" of the sphaleron process. As 
it was already mentioned, those have been found in [32] 
for COS sphalerons: they have simple exponential form. 
Lacking bound state wave functions in momentum repre- 
sentation, we can use the convolution with thermal state 



Pt 



2n t (T in ) 



f 



dp 



(2pp) 2 exp(-2pp) 
exp[(p - fi t ) /Ti n ] + 1 



(40) 



Factor 2 in denominator appears is because we only inte- 
grate over momenta directed inward, and n t (Ti n ) is the 
thermal quark density. Note that mp ~ 4 in numerical 
model wc discuss, so the corresponding effective temper- 
ature of the doorway state is l/(2p) ~ m/8 ~ 30GeV 
is softer than the actual internal Tj„ ~ 200 GeV. This 
mismatch would somehow reduce the preexponent. On 
the other hand, one can further optimize the rate, select- 
ing optimal (higher than average) momenta of the top 
quarks, to affect even the exponent. 



V. CP VIOLATION 

Finally, let us discuss the last crucial component of the 
baryogenesis, the CP violation. As it is widely known, all 
so far observed CP violation in K, B meson decays can 
be explained inside the SM, in which it comes from the 
phase of the CKM matrix. 

(The only exception to this is a very recent result from 
DO collaboration [37] on same-sign dilepton asymmetries 
which are about -0.01, a 3.2 standard deviation from 
much smaller SM prediction (-2.3 x 10~ 4 . Its origin is 
believed to be the B s <->■ B s oscillations: in a couple years 
we will know if it withstand further data samples from 
D0,CDF and LHCb. ) 

Returning to SM and its CKM matrix, it would be 
desirable to establish a specific process and prediction, 
which will follow from it for the cosmological scenario 
under consideration. Since in this scenario two (of the 
three) Sakharov's conditions - the deviation from equi- 
librium and the baryon number violation - are in a sense 
maximized, so the required magnitude of the CP viola- 
tion in this scenario need to explain the observed asym- 
metry is minimized, to something of the order of 10~ 7 or 
so. Below we discuss whether such level of CP violation 
can actually be reached inside the SM. 



A. Optimizing the one-loop CP-odd effective action 

The magnitude of the CP violation can be derived by 
standard one-loop effective action which is logdetD, with 
the Dirac operator D including Higgs and gauge fields. 



The well-known Jarlskog argument corresponds to the 
situation in which the field mass scale 



F 



W M „ » m\ 



(41) 



is larger than that of any quarks. In this case the co- 
efficient of the effective operators must be proportional 
not only to the so called Jarlskog invariant J but also to 
extremely small Jarlskog determinant 
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(m 2 b - m 2 d )(ml - m 2 s )(m 2 s - m 2 d ) ~ J 'm\m\m 2 m 2 s l 'F 12 {42) 

because when each two same-charge quarks are degen- 
erate the CP odd phase can be rotated away to zero, 
nullifying the Cp-odd effect. Since the field inside the 
bags we discuss is of the scale F ~ lOOGeV^ * \/N, this 
is indeed the case and thus we have to conclude that CP 
odd effect inside the bag is too small. 

Let us now try to optimize the effect, by looking at 
all scales F = m t ,m b ,m c ,m s , subsequently. The corre- 
sponding expressions and numerical values (for J = 10 
and actual quark masses) are 



SL C p(m t ) 



J 



l 2 2\/ 2 2\( 2 2\/ 2 2\ 
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«3*10~ 13 GW 4 (43) 



6L C p(m b ) ~ J(m 2 c - m 2 u )(m 2 - m 2 d ) * 4 * lO' 7 GeV A (U) 



6L CP (m c )~Jm 2 (m 2 -m 2 ) 



4 * 10- 7 GeV A 



5L C p(m a ) 



4 * lO^GeV 4 



(45) 



(46) 



Thus, at this level of accuracy, the scales of b and c 
give comparable effects. There are however two more 
consideration which favor the b quark scale as the optimal 
one. 

One is related with "thermal masses" in the medium of 
the quark outside bags. Weak thermal (Klimov-Weldon) 
quark mass for left handed fermions due to scattering on 
thermal W,Z is (see e.g. [35]) 



SM 2 



•rpZ 



(47) 



which for the outside temperature T ~ 50 GeV is about 
(5 GeV) 2 . This independently points to the possibility to 
go down to the scale of the b quark mass in all fcrmionic 
propagators. 

The second consideration is specific to the hybrid sce- 
nario we discuss. In this case the light quarks are pro- 
duced not from Higgs but from weak W, Z, t decays, and 
gluons from quark annihilation: and at the initial period 
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of the electroweak phase under consideration both are 
not happened yet. For this reason we have not included 
gluonic mass in the previous paragraph. For the same 
reason we expect the mechanism due to a difference in s 
and d quark scattering on the "bags" be subleading to b 
scattering or b - t transitions. 

In order to do the actual calculation one needs to cal- 
culate explicit CP violating Lagrangian and integrate it 
over the space-time volume. As a crude estimate, let us 
add the minimal power of the coupling (4Tra w ) 2 ~ 0.2 
needed for 4 CKM matrices, and "coherent space-time 
volume" V4 outside the spot in which the gauge field is at 
the scale of mb, times the estimated 5Lcp(rrib). One can 
then conclude that for V4 ~ 1 GeV~ 4 the total effect may 
reach the needed 1CT 7 . Although this volume is compara- 
ble to space between the "spots" times the lifetime, it is 
highly nontrivial that the "coherent space-time volume" 
(in which the effect maintains the same sign) can be that 
large, as the structure of the Lagrangian is complicated 
and must include the 4-dim epsilon symbol. Needless to 
say, the corresponding Lagrangian at the m& scale is not 
yet calculated (to our knowledge). 

The closest to it is the calculation in Ref.[3S] for the 
scale below m c , which produced the following operator 
of dimension-6 

8L CP (< m c ) ~ Jmfe^^ (48) 

[z^w; x w a (w;w a + w;w+) + cc ] 

Some elements of its general structure, with 4-dim ep- 
silon (providing negative P parity) and four W vertices, 
providing the Jarlskog invariant J from four CKM matri- 
ces, are generic. This structure alone implies quite com- 
plicated product of the fields, with a necessary nonzero 
electric field. 

Averaging this operator over fields found in numerical 
bosonic simulations has been performed in ref.[36 , with 
the result four orders of magnitude larger than needed. 
It is however clear that one actually cannot apply this 
Lagrangian, since this operator was derived for a com- 
pletely different field scale F « m c , which is not the 
case in the numerical simulations in question. Not sur- 
prisingly, the obtained result is larger than the observed 
one: the dimension 6 of the field would produce very 
large contribution from the "hot spots" with very strong 
gauge fields. However, as follows from Jarlskog argument 
(repeated above) at larger field there is additional strong 
suppression due to mass differences, which will kill it as 
detailed above. 



B. CP violation in top decays: the "penguins" and 
the final state interaction 

It is well known that CP violating effects in particu- 
lar meson decays are much larger than what one would 
get from the one-loop gradient expansion effective La- 
grangians like the one discussed above. Instead of high 




FIG. 5: Two interfering diagrams in top quark decays: the 
dashed line is the W and the vertical dotted line in (b) is 
the gluon. D = b,s,d,U = c,s are generic names for decay 
products, D' is the intermediate quark in the loop which is 
summed over down flavors. 

order in weak interaction, one may find the needed four 
CKM matrices in the interference terms between low- 
est order weak interaction and strong radiative correc- 
tions, the so called "penguin" diagrams. Fig(5] shows 
two generic diagrams, which can interfere and produced 
the CP-violating asymmetries from the famous complex 
phase in the CKM matrix. As seen from the figure, there 
are six generic decays for all possible values of the final 
quarks: 

t -»■ ddu, t -> ssc, t -> ddc, t -»■ bbc; t -> ssu; t -*■ bbu (49) 

The amplitude of the first one can be written as 

A(t -* ddu) = V td V* d + V tb V* b F penguin (50) 

in which we selected D' = b and introduce the so called 
penguin factor, including in particular the propagator of 
an extra gluon and the strong coupling constants. The 
interference term between these two terms would provide 
the needed four CKM matrices, while the "time arrow" 
can be provided by a specific direction of the flow, e.g. 
leakage of bound tops into the unbound lighter quarks 
leaving the bag. (Note again, that in equilibrium there 
would be inverse processes which will cancel the asym- 
metry, but those are absent in the scenario we discuss: 
the light quarks are very much out of equilibrium and 
their density is negligible at the time we consider.) 

The ratio of the two terms generates CP asymmetry, 
which for this decay is as large as 

„ p VtbV u bFpenguin p , 

V td V ud 

Other modes have smaller asymmetry by powers of A, 
but their decay rates are larger. 

To get the average effect one has to evaluate the ef- 
fective CP-odd Lagrangian or the imaginary part of the 
product of the two terms, their interference. For all 6 
decays mentioned above 

— — ~ a w Im[VtdV ud V tb V ub ]F pengum 

~ al V X 6 F penguin « 10- 8 (52) 

where we have used the Wolfenstein parameterization (to 
A 4 accuracy) of the CKM matrix, and include "penguin 
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suppression" factor ~ 1/3. The net effect of this La- 
grangian would be a difference between the number of 
bound tops and antitops in the bag, and with "top recy- 
cling" leading to asymmetric sphaleron transitions and 
baryon asymmetry. 

(Note that the quark mass differences entering the Jarl- 
skog determinant argument of the previous subsection, ft 
resigns in the summation over the D' quarks. If all three 
Penguin diagrams would be exactly the same, one would 
get the following combination 

V tb V; ib + V ts V: s + V td V: d (53) 

which is zero because of the orthogonality condition of 
the two rows. Here however where different quark masses 
and "external conditions" (in form of the external ther- 
mal mass in the propagators of the D' quarks) should 
help to produce different coefficients. Since the scale of 
such external effects are of the order of F ~ few GeV, 
they clearly separating the b penguin term from two oth- 
ers. As in the previous section, the remaining suppression 
of the order of (m s - m d )/ F may yet appear. ) 

Most CP in meson decays are due to mixing of the 
neutral meson-antimeson states: but this is not needed. 
Top/antitop decays are like those of the charged meson, 
which will have a CP effect due to the final state interac- 
tion, or "strong phases" as they are usually called. The 
first and second diagrams should be "dressed" by the 
gluon exchanges between the quarks (not shown in Fig.), 
it will result in some phases 5 dlrectl S penguin ~ 0(a s ), 
which are the same for quarks and antiquarks. They 
are not the same for different final states, for example if 
the quark pair produced from a gluon in a penguin has 
small invariant mass, it corresponds to a pair produced at 
relatively large distances from the original quark, while 
the "direct" decay diagram is near-pointlike because of 
large W mass. Standard arguments show that there is a 
difference in the decay rates of particle and antiparticle 

\ A (t -* f)\ 2 ~\ A (t -> I)? ~ sin(8 direct - S penguin )sin(6 C p) 

(54) 

Only recently the first such asymmetry has been observed 
in a difference between the width T(B + -> poK + ) decay 
and that of its CP conjugate. So, there is no doubt that 
partial widths of the t and t decays can differ. While the 
CPT theorem forbids the difference in their total width, 
for time dependent bags well out of equilibrium one may 
still hope that it would be possible to go around strict 
CPT it and generate a difference in t and i population, 
perhaps even of the magnitude comparable to the ob- 
served asymmetries in its decay modes. 



C. Baryon charge separation during the bags 
formation 

Another possible manifestation of the CP-odd phe- 
nomena is possible separation of the quarks vs antiquarks 



already during the process of bag formation, in the ini- 
tial very off-equilibrium stage. If it happens, then the 
amount of tops and antitops in the bag would already be 
different, leading to an asymmetry in sphaleron transi- 
tions which "recycle" those tops, as detailed above. 

Farrar and Shaposhnikov [38] were the first who dis- 
cussed it in detail, in a specific case of strange quarks and 
near-equilibrium bubbles at T = Tc, assumed at the time 
to be the first order transition. Although we would use 
this phenomenon in a quite different setting, we would 
still call it "the Farrar-Shaposhnikov (FS) mechanism". 
Its essence is the interference between the (CP-even) 
scattering/reflection on the bag wall and the (CP-odd) 
transition via the CKM matrix into a different flavor, re- 
flection, and return to the original flavor. Let us remind 
the reader how it works. Let us start with flavor f\. if it 
is reflected from the bag wall it gets the phase 6\ . Alter- 
natively it can turn into the flavor 2 and scatter with the 
phase 62 , picking on the way CP-odd phase from complex 
CKM matrix elements. The total baryonic current is a 
difference between those for quarks and antiquarks 

IgWi + (j e i8 2 +i8 C p | 2 _ LiSi + (j e iS 2 -iS C p | 2 

~ sin(6i - 52)sin(6cp) (55) 

The essential part of the argument is that two flavors 
under considerations should have different scattering on 
the bag, as the effect disappear at 5i = 82- This condi- 
tion is especially difficult for light quarks, as they only 
weakly interact with the Higgs and thus the bag. FS 
has argued that the effect should be nonzero for s rela- 
tive to d quarks, but only in a narrow strip of energies 
near the level crossings, where s quarks are all reflected 
(large 8\) while d are transmitted (negligible S 2 )- As a 
result, their estimate for the FS effect is proportional to 
a small parameter corresponding to the relative width of 
this region, (m s - m d )/T ~ 1CT 3 . 

Another important part of their paper (which we would 
not remind here) is discussion of the unitarity and CPT 
constraints on possible effects, leading to a standard con- 
clusion that in equilibrium (e.g. for a bubble wall at rest 
with the medium) the effect vanishes, as it indeed should 
from general considerations [T]. 

Let us now point out large differences which exist be- 
tween the matter properties discussed by FS and that 
in the hybrid scenario under consideration. First of all, 
there are drastic differences in its chemical composition. 
Equilibrium matter at electroweak scale has about hun- 
dred effective degrees of freedom, including all quarks and 
leptons, as well as gluons. Therefore effective fermionic 
thermal masses are dominated by the gluonic part, with 
a scale ~ T/2 ~ 50 GeV. In the hybrid scenario there was 
not yet enough time to generate light quarks, leptons and 
gluons: the only quanta produced are weak gauge bosons 
W, Z, Higgs and the top quarks, all due to their large cou- 
pling to (strongly fluctuating) Higgs field. The spatial 
distribution of tops is very space and time-dependent, as 
they flee the regions with large Higgs VEV and are col- 
lected into the bags we discuss. While doing so, they can 
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also migrate, via the CKM matrix, into all other flavors, 
driving the FS mechanism and producing different top 
and antitop population in the bags. 

For simplicity, let us imagine the bag to be very large 
and its boundary flat, separating the space into two halfs, 
with symmetric (v « 0) and broken (v i» Vq) phases. The 
flow, into symmetric phase, is driven by the difference in 
the chemical potential. As we have discussed throughout 
this paper, the temperature is also not constant over the 
divide, with T ~ 50 GeV in the broken phase outside the 
bags. Tops in the bags (symmetric half) with energy less 
than the top mass outside E < M t (T out ) cannot leave the 
bag and are reflected back if collide with the boundary: 
thus large St- If they however migrate into the c, u quarks, 
those hardly see a bag at all, thus 5 C ,5 U « 0. Additional 
small parameter appearing in FS estimate is the part of 
the energies in which c, u contributions do not cancel is 
(m c -m u )/T ~ 1/50, not so small. Together with Jarlskog 
J from CKM, it will generate asymmetries of the order 
of 10~ 6 or so, which is right in the needed ballpark. 

(A word of warning: The original FS estimate of the 
effect has been criticized in refs [55] . One of the points 
is that additional suppression appears from the widths 
of quarks coming from in media scattering, and we are 
not sure what is the final word in the polemics about 
it. Anyway, this criticism is much less relevant to the 
scenario discussed, as larger scale m c is substituting m s 
in FS, which is rather close to - much reduced relative 
to equilibrium matter - rescattering widths.) 



VI. SUMMARY AND DISCUSSION 

In summary, we have considered the role of the W - Z- 
top bags in baryogenesis. We have used, as a specific 
example, numerical simulations for the hybrid preheating 
scenario from Ref . [29] . 

(i) We first reproduced the main results of those simu- 
lation using the W-Z bags and COS sphaleron, including 
qualitatively the sphaleron size and rate found numeri- 
cally. 

(ii) We then discussed the role of the top quarks and 
found them to be quite important. They are produced 
via HH -* it more effectively than W. They are collected 
into the bags together with the gauge quanta, presumably 
producing bags of larger size, with improved lifetime. We 
however found that the mechanical stabilization condi- 
tion (which makes bags with tops cooler than pure W- 
bags) basically makes the sphaleron exponent constant, 
independently of how many top quarks are in the rate. 

(iii) In addition, there is an interesting effect of "top 
quark recycling" . This effect helps the transition, pro- 
viding significant (~ 1/3) fraction of the energy needed 
to reach the sphaleron mass. The recycling of 3 tops 
increases the rates by another factor 20. 

(iv) We then provide arguments that optimal CP vi- 
olation should happen either due to high order weak 
interaction in the area between the bags, at a scale of 



nib ~ 5 GeV, or in top decays, or in top selection during 
bag formation due to Farrar-Shaposhnikov effect. Al- 
though we are not able to provide specific calculation 
of any of the three possibilities at this time, crude esti- 
mates show that any one of them can potentially provide 
a magnitude of 10 which, together with the estimated 
sphaleron rate, would generate the baryon asymmetry of 
the observed magnitude from known CKM-base mecha- 
nism alone. (If the DO result on much larger CP viola- 
tion in semileptonic B s decays would stand in time, one 
should of course also look at its origins as well.) 

Let us finish this paper with a historical analogy. We 
do see stars on the sky which shine: yet basic nuclear 
physics was not enough to explain why they can do 
so, falling short by many orders of magnitude. It took 
decades of patient and specific work to uncovered the 
nontrivial ways in which the star cycles work. Those cru- 
cially depend on some catalysts, whose presence is very 
small but crucial. Perhaps the W - Z-top bags are also a 
kind of a catalyst, or a doorway state, helping to increase 
the probability of the baryon number violation and of the 
CP violation. 
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Appendix A: The coupling constants and parameters 

Let us remind the reader the hierarchy of the coupling 
constants, as it reflects the order in which various pro- 
cesses should be considered. 

The largest coupling constant is still that of strong 
interaction, 

2 

— = a a (100 GeV) * 0.11 (Al) 

47T 

but all the original bosonic fields of the electroweak the- 
ory W, Z, H can only interact with strong interaction via 
quarks, so it has to play so to say the secondary role. 

Quite close to it is the Yukawa coupling of the top 
quark to Higgs, if written in the same form it is 

^ = a t * 0.08 (A2) 

4-7T 

and this would be responsible for most of the effects 
we will discuss below. Other Yukawa couplings are 
much smaller and can be neglected: for example for the 
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next b quark such combination is smaller by the factor 
(m b /m t ) 2 m 7.4 * 1CT 4 . 

Weak interaction coupling is naturally smaller 



.03 



(A3) 



in the standard model: note however that in numerical 
simulation we discuss a much smaller value has been used, 
presumably for some methodical reasons. We will not in- 
clude weak interactions except in top quark decays. Like 
in many other works, we will ignore electromagnetism 
(the U(l) gauge part), putting the Weinberg angle to 
zero: thus mjy = mz- 

It is well known, the Higgs mass (or self-coupling A) 
remains the last unknown parameter of the SM. Our pre- 
vious studies of the top-balls put the Higgs mass to a 
round number Mh = 100 GeV. 

However, since we try to explain numerical bosonic 
simulations [29] , let us mention that their choices for the 
Higgs mass and VEV are different. Details can be seen 
in their original work, for example the particular set of 
parameters, called Al in |29| for which the pictures we 
used are based have rather heavy Higgs mass 



A = g 2 /2 = 0.00675, m H /m w = 4.65 



(A4) 



and also a different the value of the electroweak coupling 
is g w = 1/20. The Higgs VEV is still v = 246 GeV as 
in real world. The mass parameter of the Lagrangian 
related to physical Higgs mass by to = m H / v2 is chosen 
as a basic mass/energy/length/time unit, which we will 
also adopt in what follows. With realistic Mw it makes 
77i = 264 GeV, not far from v. 



Appendix B: Strong annihilation 

The differential cross sections for these two QCD re- 
actions are well known, here we present them ignoring 
quark masses |H] : 



datt^qq _ ^a 2 s 



dt 



[1 + u 2 /s 2 ] 



da 



tt-*gg 



4ira 2 ,t 2 + u 2 



dt 



Integrated it we get[ 



a tt^qq ~ a tt- 



ut 



39 



)( 



8 6ut 



0.022 



(Bl) 



(B2) 



(B3) 



Thus the contribution of those two processes to the decay 
rate of the top quark can thus be written as 



2 * 0.022- 



/(Pl)/(P2) 

(2tt) 3 (2tt) 3 2pip 2 [l-co»(e 12 )] 



/ <&/(Pl) 



(B4) 



where f(p) is the thermal distribution and 9i2 is the 
angle between momenta of two colliding tops. Assum- 
ing that logarithmic angular integral is regulated by 
(1 - cos(8 min )) « 0.1 we get the following simple answer 
for the strong annihilation time 



7* = 1/Tft « 200/T 



(B5) 
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